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Abstract 

The Wilkinson Microwave Anisotropy Probe (WMAP) experiment has detected reionization at 
the 5.5a level and has reported a mean optical depth of 0.088 ± 0.015. A powerful probe of 
reionization is the large-angle EE polarization power spectrum, which is now (since the first five 
years of data from WMAP) cosmic variance limited for 2 < Z < 6. Here we consider partial 
reionization caused by WIMP dark matter annihilation, and calculate the expected polarization 
power spectrum. We compare the dark matter models with a standard 2-step reionization theory, 
and examine whether the models may be distinguished using current, and future CMB observations. 
We consider dark matter annihilation at intermediate redshifts (z < 60) due to halos, as well as 
annihilation at higher redshifts due to free particles. In order to study the effect of high redshift dark 
matter annihilation on CMB power spectra, it is essential to include the contribution of residual 
electrons (left over from recombination) to the ionization history. Dark matter halos at redshifts 
z < 60 influence the low multipoles I < 20 in the EE power spectrum, while the annihilation of 
free particle dark matter at high redshifts z > 100 mainly affects multipoles I > 10. 
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I. INTRODUCTION 



Studies of quasar spectra indicate that the Universe is highly ionized up to a redshift 
z ~ 6. In the standard approach, early stars, early galaxies, quasars, etc. are thought to be 
the sources responsible for reionizing the Universe. The free electrons that exist after reion- 
ization scatter microwave photons, resulting in an increase in the degree of polarization of 
the cosmic microwave background (CMB) radiation. The Wilkinson Microwave Anisotropy 
Probe (WMAP) has confirmed the reionization of the Universe and has computed the optical 
depth due to scattering of CMB photons with free electrons [IH5]. 

There is considerable uncertainty regarding the extent to which early stars reionize the 
Universe. The epoch of formation of the first stars, their masses, their abundance, etc are 
model dependent factors that influence the ionization efficiency There are also questions 
regarding the nature of the first stars. In some models with weakly interacting massive 
particles (WIMPs) being the dark matter, the earliest generation of stars were dark stars 
[6] characterized by a low surface temperature, and very high mass. In these theories, it is 
unclear how the next generation of stars formed, and whether they were able to partially 
reionize the Universe. Here, we will consider dark matter annihilation to be a source of 
reionization. 

Several authors have studied the impact of particle decay and annihilation on reioniza- 
tion j7]. Annihilating dark matter with mass m x ~ 10 — 100 GeV may be detectable by 
future CMB polarization measurements [E] . More recently, there has been an interest in con- 
straining the properties of dark matter models using the WMAP measured optical depth. 
In previous work [HI [TO], we studied dark matter halos with low WIMP masses, with the 
restriction that they be consistent with the measured optical depth. Heavier dark matter 
particles m x > 100 GeV in halos could reionize the Universe for favorable particle physics 
parameters [11], and could possibly account for the positron excess in high energy cosmic 
rays found by ATIC [12] and PAMELA [13] . Reionization by annihilating dark matter for a 
variety of particle and cosmological parameters was presented in [H] ; reference [15] studies 
the allowed regions in the dark matter mass-cross section plane. More recently, [TB] have 
analyzed a combination of recent CMB data sets for evidence of dark matter annihilation. In 
this article, we investigate whether it is possible to distinguish standard reionization theories 
from more exotic dark matter scenarios, using CMB observations. 

Scattering of free electrons by CMB photons causes additional polarization at scales ~ the 
horizon at reionization, resulting in excess power at low multipoles compared to a Universe 
that is not reionized. Let us denote the WMAP measured optical depth as r. We note that 
this is generally different from t(z), the total optical depth up to redshift z. This is because 
the WMAP experiment may not be sensitive to large z, and also because of the effect of 
residual electrons included in r(z). The large angle EE polarization oc t(z) 2 is an excellent 
probe of reionization, but is seen at high significance by WMAP only for Z < 10 [31 IS] ■ With 

1 year of data, the WMAP experiment detected a reionization signal in the TE polarization 
power spectrum, implying an optical depth r = 0.17 ± 0.04 pQ. This value was shown to be 
too large by the WMAP 3- year analysis [2]. With 3 years of data, WMAP obtained a value 
for r = 0.10 ± 0.03 using EE data alone, and r = 0.09 ± 0.03 using EE, TE, and TT data 
PQ . The 3-year measurement of r has been confirmed by the 5-year analysis [31 H] • The 5- 
year WMAP measured EE power spectrum is largely cosmic variance limited at multipoles 

2 < Z < 6 [3J. The WMAP 5-year analysis provides a mean value r = 0.087 ± 0.017 
using all WMAP data [I]. With BAO and SN data included, the measured optical depth 
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is t = 0.084 ± 0.016 [T7]. Assuming an instantaneous reionization model, WMAP reports 
a reionization epoch z* = 11 ± 1.4 using WMAP data alone, and = 10.8 ± 1.4 using 
WMAP+BAO+SN data |17] . These results are in excellent agreement with the recent 7- 
year result r = 0.088 ± 0.015 for sudden reionization and r = 0.087 ± 0.015 for a fit which 
allowed a varying width of a one step reionization scenario. The low I polarization data 
(EE) alone, implies a detection of reionization at the 5.5a level pjBJ. 

In Section II, we provide a brief discussion of dark matter annihilation and the inverse 
compton scattering process. We compute the energy absorbed by gas at a given redshift, 
and the expected ionization and temperature histories. In Section III, we consider dark mat- 
ter models with different particle masses and concentration parameters, and calculate the 
expected polarization power spectra using the camb program [19] with appropriate modifica- 
tions. We first consider dark matter annihilation by free (gravitationally unbound) particles 
at high redshifts z > 100 with and without residual electrons. Ionization in this regime 
directly probes the particle physics properties of the model. We then consider ionization 
at all redshifts, including WIMP annihilation in halos, and compare with a standard 2-step 
reionization model. Finally, we present our conclusions. We used the parameters Q = 1, 
Q h h 2 = 0.022, Q dm h 2 = 0.12, h = 0.7, A s = 2.3 x 10~ 9 , n s = 1, w = -1, T cmb = 2.726 K to 
make our plots. 



II. PARTICLE ANNIHILATION IN HALOS. 

The number of dark matter annihilations per unit volume and per unit time, at a redshift 
z is given by: 

dN ann (a a v) 2 

W (Z)= 2K P * (4 () 
(a a v ) is the annihilation cross section times the relative velocity of the particles, averaged 
over the WIMP velocity distribution. m x is the particle mass, and p x (z) is the dark matter 
density at z. We have assumed that two WIMPs participate in an annihilation process 
producing particles with total energy 2m x c 2 per annihilation. Before the epoch of halo 
formation, p x takes the simple form: 

Pl free (z) = (i + zf p 2 c n 2 dm , (2) 

where p c is the critical density, and f2dm is the dark matter fraction today. 

As the annihilation rate is proportional to the square of the WIMP density, it is enhanced 
by the formation of structure in the Universe. Hierarchical structure formation starts at a 
redshift z ~ 60 [201 EE]- The boost in density due to clustering of WIMPs is partially 
countered by the small halo volume. Nevertheless, p x .haio eventually becomes much larger 
than p x ,free (unless the concentration parameter is very small). p x ,haio is given by: 



P Xlh aloW = (l+^) 3 / dM d ^{M lZ ) 



r 20Q 

dr4nr p h (r) 



o 



(M,z). (3) 



^haio is the comoving number density of halos. The volume integral over the halo is a function 
of both halo mass M and redshift z. Here we have assumed that halos are truncated at r2oo, 
the radius at which the mean density enclosed equals 200 times the background density at 
the time of halo formation, ph(^) is the halo density at a distance r from the halo center. 
Eq. ([3]) ignores halo-halo interactions, as well as interactions between free dark matter and 
dark matter in halos. 
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A. Ionization and heating. 

WIMP particle annihilation generally results in the release of photons, charged particles, 
and neutrinos. Let dN 7 /dE 7 (z) be the number of photons produced per unit energy per 
annihilation at z. dN^/dE^ includes both prompt photons released by WIMP annihilation, 
as well as lower energy photons produced as a result of processes such as inverse compton 
scattering (ICS henceforth), bremsstrahlung, etc. The energy absorbed at a redshift z per 
gas atom £(z) can be calculated as [TO] : 

(4) 

In Eq. Q, we expressed the time variable dt in terms of dz using the relation dz = 
—dt(l + z)H(z), where H(z) is the Hubble parameter at redshift z. The cubic term accounts 
for the expansion of the Universe in the time it takes for the photons to travel from z' to 
z. E' = Ej(l + z)/(l + z'), which accounts for the redshifting of photon energy. We have 
assumed that the photons are produced at the redshift of WIMP annihilation z' . This is 
exact for prompt photons, and a good approximation for ICS photons produced by charged 
particles interacting with the CMB. E\ and E 2 are suitably redshifted energies. The last 
term in square brackets accounts for scattering with gas atoms at z. c is the speed of light, 
while a is the cross section. The exponential term accounts for the scattering of photons as 
they travel from z' to z. k(Ej) is the optical depth for photons of energy E J} given by the 
expression 

- z „ )H ^ cn(z>(£ 7 ). (5) 

Let r]i on (z) and rfa ea t(z) be the fractions of energy that go into ionization and heat- 
ing respectively. These fractions have been calculated by [221425] . For simplicity, we set 
^ion = ??heat = 1/3, with the remaining contributing to collisional excitations and the Ly-a 
background. The ionized fraction x- lon (z) and the gas temperature T(z) may be calculated 
by solving together, the two equations: 

- (1 + z)H(z) ^ Xl ° n Z = fi [1 - x ion (z)} i] ion (z)£(z) - n(z)x 2 lon (z)a(z) 

-(l + z )H(z)^& = -2T{z)H{z) + 2r] XT iZ) Z(z) + XUZ) [T f | - T(Z)] • (6) 

dz 3/bb tc[z) 

fi ~ 0.07 eV -1 is the inverse of the average ionization energy per atom, assuming 76% H and 
24% He, neglecting double ionization of Helium a is the recombination coefficient, T 7 is 
the CMB temperature, and t c is the Compton cooling time scale ~ 1.44 Myr [30/(1 + z)] 4 . 
The last term in the temperature evolution equation accounts for the transfer of energy 
between free electrons and the CMB by compton scattering EZ]- We used x ion 1 and 
ignored the Helium number fraction in the temperature coupling term. 

B. Inverse compton scattering with the CMB. 

As mentioned previously, WIMP annihilation results in the production of charged parti- 
cles, photons, and neutrinos. The charged particles rapidly lose energy by inverse compton 
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scattering with CMB photons, producing a spectrum of lower energy photons of energy E ics . 
These ICS photons are far more efficient in ionizing and heating the gas than the prompt 
high energy photons [Uj. The physics of inverse compton scattering has been reviewed by 
[25] and more recently by |29j . 

At high redshifts, the energy loss is extremely rapid, with a ~ 1 GeV electron losing 99% 
of its energy in < 0.02 Myr at z = 50, and ;$ 2 years at z — 500. The average energy of the 
ICS photons is given by 



(£ ics > = - 7 2 <6> 



l.eMev(§§) 



l + z 
501 



(7) 



where E- lcs is the energy of the ICS photons, 7 is the Lorentz factor of the electron, and e is 
the CMB energy before scattering. We used the mean value (e) « 2.7fc b T 7>0 (l + z) in Eq. |7) 
assuming a black body spectrum. 7 is time dependent since the electron is losing energy. 
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FIG. 1: Inverse compton scattering of a single electron of energy 1 GeV with the CMB at z = 500. 
The curves show the photon spectrum at times when 1%, 50%, 90%, and 99% of the initial electron 
energy has been lost to the CMB. 

The ICS photon spectrum of a 1 GeV electron interacting with the CMB at z = 500 is 
shown in Fig. [Tj The four curves in Fig. [T]are drawn for 1%, 50%, 90%, and 99% energy 
loss. Initially, the ICS photons are peaked at ~ 1 MeV, in accordance with Eq. ([7]). After 
90% of the initial energy is deposited into ICS photons, further scattering leads to new ICS 
photons with average energy ~ (0.1) 2 MeV. These photons add to the high energy ICS 
photons produced at earlier times. The result is a broad spectrum with plenty of ionizing 
photons. 

At high redshifts z ~ 500, the few ~ MeV energy photons may be expected to scatter with 
gas atoms, and produce charged particles by pair production and compton scattering. These 
high energy particles will then lose energy to the CMB, adding to the ICS photon spectrum. 
Thus, we expect that at high redshifts, some of the high energy photons are "recycled" into 
lower energy bins. The authors of [25] find that at low redshifts 10 < z < 50, a sizeable 
fraction of the electron energy is converted to X-ray photons, which may free stream to the 
observer. The lower efficiency of inverse compton scattering at low redshifts must be taken 
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into account when accurate results are required. In this article, we will be focused mainly on 
distinguishing different scenarios using CMB data, and will assume that ICS photons are of 
sufficiently low energy to effect ionization and heating. The detailed calculation of photon 
energies for arbitrary initial electron spectra is left to future work. 

Assuming ionization and heating are primarily due to ICS photons with E ics < MeV, we 
may make the assumption cr(E 7 ) « <tt, where <tt is the Thomson cross section, k is then 
energy independent, and the energy integral is simply 

where f em is the fraction of energy in the form of photons/charged particles, i.e. with elec- 
tromagnetic interactions, and the energy per annihilation = 2m x c 2 . 



III. OPTICAL DEPTH, REIONIZATION, AND CMB POLARIZATION. 

The optical depth due to scattering of microwave photons with free electrons is given by 
the expression 

r( Zl ,z 2 ) = ca T ^ H{z){1 + z) *U*M*) 0) 

~ 0.04 for Z\ = 0,z 2 = 6, assuming x ion = 1, doubly ionized Helium for z < 3, and H(z) = 

Ho\J^m(^ + z) 3 + Qa- This value of r(0, 6) is less than one half of the WMAP measured 
optical depth. Thus, the Universe was at least partially ionized at redshifts z > 6. As 
mentioned in the Introduction, the WMAP measured EE power spectrum is most sensitive 
to the low / multipoles. The EE power spectrum data at multipoles I < 10 is sufficient when 
reionization is due to astrophysical objects at redshifts z ^ 25 jl]. In the standard theory, 
this is the case. However, it is not the case when dark matter annihilation is the source of 
ionization. We now consider partial ionization by dark matter annihilation at intermediate, 
and high redshifts. 



A. Ionization at very high redshifts z > 100: Annihilating free dark matter 

At redshifts z > 100, Eq. Q may be further simplified since photons are absorbed close 
to the redshift at which they are released. We can approximate the exponential e~ K in Eq. 
Q as a delta function normalized so that A/" / dz'e~ K = 1: 

co^n(z) 

which makes £ independent of the cross section a. We have verified by performing the full 



calculation that £ is indeed independent of a when a is close to ctt- Using Eq. (10) and Eq. 
Q, we obtain the simplified formula, valid at high redshifts z > 100: 



cm 



m x n(z) 



. X in-3^L + (f*A M/(3xlQ- 26 cmVs) 

Myr I 501 J \2/3J m x /(100GeV) 1 ' 
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Eq. (11) shows that dark matter ionization at redshifts z > 100 can directly probe the 
quantity / em (fj a ?;) A%> without any of the complications of the halo model. / em ~ 2/3 
assuming two thirds of the energy goes into charged particles and photons. To see the 



physical significance of Eq. (11), let us multiply the result with a typical time scale = 
l/H{z) sa 2.4 Myr [501/(1 + z)f /2 . Assuming 1/3 of this energy is available for ionization 
and heating respectively, we get (£/3)£h ~4x 10 -3 eV at z = 500. Since the mean energy 
of the CMB at this redshift = 2.7&bT 7i o(l + z) ~ 0.3 eV, we might not expect dark matter 
annihilation to affect the gas or CMB temperature significantly. Assuming the same energy 
is available for ionization, and an ionization potential of 13.6 eV, we may expect a fraction 
4 x 10~ 3 /13.6 ~ 3 x 10 -4 (100 GeV/m x ) of the atoms to be ionized (without recombinations). 
Since this number is comparable to the residual electron fraction ~ few xlO -4 , we expect 
dark matter annihilation to contribute significantly to ionization, particularly for small m x . 
We also expect the increased ionized fraction to better couple the gas and CMB temperatures 
by compton scattering. 

Let us consider 4 different dark matter models with masses m x = 10, 50, 100, and 500 
GeV. Let x- m i t be the value of x ion at our starting redshift z = 1000. If a; init is too small, 
the first term in Eq. ^ will drive x ion up to the appropriate value. Similarly, if x init is 
too large, x ion decreases due to the second term in Eq. This ensures that the value of 
x; on far from the initial redshift z = 1000 is independent of the chosen value of Xmit- Let us 
model Xion as a power law in (1 + z), valid for z > 100, i.e. 

/ 1 ~f~ z \ n 

X ion (z) = Xion(^) ( : ; ~ ) , (12) 



1 + -2* 

where is some reference redshift far from z = 1000 where the 3 curves are nearly coincident. 
The exponent n is determined using the value of Xi on at another redshift z where the 3 curves 



are close to one another. We may then extend Xi on using Eq. (12) to obtain the "power law 
fit" initial value Xj n i t . 

In Fig. [2j the solid(red) line denotes the power law fit to the 3 curves shown by dashed 
lines. For small WIMP masses, the curves become coincident quickly, while for larger masses, 
they approach each other only gradually. Once the value of n and [z*, Xi 0n (z*)] are known 
for any given m x , Eq. (12) may be used to compute r(60, 1000), the contribution to the 



optical depth for 60 < z < 1000. 



B. Ionization at very high redshifts z > 100: Residual electrons 

The above discussion neglects contributions from residual electrons from the formation 
of helium and hydrogen atoms, which continues well after photons decouple (as there are 
many more photons than baryons/electrons in the Universe). Due to residual free electrons, 
Xion never falls below ~ 10 -4 — even in the absence of ionizing sources before z ~ 10. 

When residual electrons are included, the importance of the contribution to the optical 
depth made by dark matter annihilation is considerably reduced because the recombination 
term in Eq. ^ oc xf on is dominant at high redshifts. Let tr be the contribution to the optical 
depth made by residual electrons. To estimate tr, we use the recfast code described in 
[27] . and implemented in the camb program [19]. The contribution of residual electrons to 
the optical depth has been studied by [32J , who find that residual electrons may contribute 
~ 0.06 to the optical depth between redshifts 7 < z < 700. The optical depth contribution 
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FIG. 2: Estimating the initial value of Xi on for the different dark matter models in the absence of 
residual electrons. The solid line is the power law fit, while the dashed lines are drawn for different 
values of x m \t ■ 

tr(z) increases sharply for z > 700. Table 1 lists the values of r(60, 1000) obtained using 
the power law fit of x ion for WIMP masses m x = 10, 50, 100, and 500 GeV respectively, 
with and without residual electrons. As seen from Table 1, the effect of residual electrons 
is to reduce the impact of the optical depth contribution of dark matter annihilation to 
low values except for the case of very light particle masses, i.e. dark matter annihilation is 
ineffective when the residual electron fraction ^> £(z)/[(13.6 eV)H(z)]. Nevertheless, dark 
matter annihilation always increases the total optical depth. Thus, one needs to carefully 
consider the recombination history when placing constraints on dark matter annihilation at 
high redshifts. 

C. Reionization and CMB power spectra. 

We now consider dark matter models at both high and intermediate redshifts, and com- 
pare these with a standard, two step baryonic reionization model (inspired by simulations 
as e.g. [23]), referred to as Model A. We use the term 'baryonic' simply to indicate that 
reionization is entirely due to luminous astrophysical sources, as opposed to dark matter 
annihilation. We consider a 2-step reionization scenario for Model A. From z = 25 to z = 6, 
we assume a linear increase of x- lon (z) with z, with full ionization for z < 6. Xi on (z) is 
smoothly patched on to the value of Xi on for z < 6 and z > 25 using a fitting function of the 
hyperbolic tangent type. The particular form of Xi on , as well as the fitting function for this 
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Values of r(60, 1000) for different dark matter masses. 



m x (GeV) 


10 


50 


100 


500 


r(60, 1000) (r R = 0) 


0.070 


0.033 


0.024 


0.011 


r(60, 1000) - r R 


0.028 


0.007 


0.004 


0.001 


r(60, 1000) 


0.429 


0.407 


0.404 


0.401 


n 


0.295 


0.253 


0.234 


0.171 


ZionO*) x 10 4 


8.9 


4.2 


3.1 


1.5 



TABLE I: Estimated values of the optical depth for the different dark matter models with and 
without residual electrons. The row r(60, 1000) (tr = 0) is due to dark matter annihilation only, and 
assumes zero residual electron contribution. The next two rows include residual electrons computed 
using the recfast code. r(60, 1000) — tr is the contribution due to dark matter in the presence 
of residual electrons, while t(60, 1000) is the full optical depth, which includes contributions due 
to both dark matter and residual electrons. Most of the contribution to r(60, 1000) comes from 
z > 700. n and Xi on (-z*) are defined in the text, « 400. The WMAP measured value of optical 
depth is mostly due to ionization at late times [z < 25) and does not include the contribution due 
to residual electrons. 



model are chosen primarily for convenience, and are not motivated by any particular theory 
of star formation. 

We also consider 3 models motivated by dark matter annihilation. In these models, 
reionization is due to both dark matter and luminous sources. For z < 25, we assume the 
same reionization scenario considered in Model A. For intermediate redshifts 25 < z < 60, 
partial reionization is due to dark matter annihilation in halos. For z > 60, free dark matter 
particles play the dominant role. In the intermediate redshift range 25 < z < 60, where 
dark matter halos are the primary source of ionization, we assume a Navarro-Frenk- White 
[54"] type density profile, and a Sheth-Tormen [32] type distribution of halo masses. The 
concentration parameter of the halos is treated as a free parameter. DM Model $T has a 
particle mass m x = 100 GeV, and a concentration parameter c = 15. DM Model #2 has m x 
= 10 GeV, c = 5, and DM Model #3 has m x = 10 GeV, c = 15. The optical depth at low 
redshifts r(0, 25) is set to a value of 0.087 for all models, consistent with the measurements 
reported by WMAP. 

To quantify the difference between two models X and Y, let us compute the a{l) deviation 
for the E mode polarization, denned as: 




°V) = ^7^- ■ (13) 



Ci = Cf E is the 2 point EE polarization power spectrum for multipole I. AC/ = AC EE is 
the la cosmic variance error at I: 

< AC ') 2 = WTi c ' • < 14 > 



which neglects instrument noise and assumes 100% sky coverage. In general, Cf in Eq. (14) 
should be replaced by 

r) 2 = \\cf A cr+(cr) 2 



(15) 
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FIG. 3: Comparing the 3 dark matter models with baryonic Model A. For z < 25, the dark 
matter models are identical to Model A, and have an optical depth r(0, 25) = 0.087. For z > 25, 
the dark matter models feature partial ionization due to halos and free particles. Also shown for 
comparison is a 'sudden' reionization model, (a) shows the ionization history for the different 
models, (b) shows the respective polarization power spectra (EE), with la cosmic variance error 
bars computed using the power spectrum of Model A. (c) shows the a (I) deviations of the different 
models (EE power spectra), all compared to Model A. The solid (black) curve denotes the 'sudden' 
reionization model, while the 3 dashed curves are drawn for the 3 dark matter models. DM Model 
#1 is indistinguishable from Model A. 



where A and B could stand for T or E. Instrument noise is not negligible in real experiments. 
For example, instrument noise for the Planck experiment is expected to become comparable 
to cosmic variance at I = 20 |30j . However, in this article, we use Eq. (14) and consider a 
cosmic variance limited observation. 

Fig. |3^a) shows the ionization histories of the different models. The baryonic model A 
is shown by the solid (red) curve. For comparison, we also include a 'sudden' reionization 
model, shown by the solid (black) curve. In the sudden reionization model, the Universe is 
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ionized below z ~ 11. The three dark matter models are shown using dashed lines. We note 
that the two models with m x = 10 GeV (DM Model #2 and DM Model #3) have larger 
X; on at high redshifts, compared to the model with m x = 100 GeV (DM Model #1). Also 
of note is the effect of the concentration parameter c. (b) shows the respective polarization 



power spectra. The error bars are cosmic variance only, and are calculated using Eq (14) 
for model A. The dashed curves are drawn for the 3 dark matter models, (c) shows the a(l) 
variation (comparing EE power spectra) between the 3 dark matter models and Model A. 
The solid (black) curve denotes the sudden reionization model compared with Model A. The 
difference between DM Model #1 and Model A is nearly zero. DM Models #2 and #3 show 
a larger deviation and may be distinguished from Model A at high confidence with cosmic 
variance limited data up to Z = 20. With cosmic variance limited data up to I = 10, the 
dark matter models are not easily distinguishable from the standard reionization scenario. 
Considering that WMAP polarization data is cosmic variance limited only for 2 < I < 6, and 
the signal is not significant for I > 10, current large angle polarization observations cannot 
distinguish between different reionization histories, but future experiments may be able to 
do better. [31] have considered separately, ionization at redshifts 6 < z < 10 and z > 10 
and have constructed exclusion contours given the WMAP data, as well as predictions for 
Planck and CMBPol. However, the fiducial ionization history considered in [31] does not 
include ionization by dark matter at redshifts z > 100. We conclude that in the absence of 
polarization data for I > 10, it is hard to constrain dark matter annihilation (m x > 10 GeV, 
thermal relic cross section), or differentiate between different reionization models using large 
angle CMB polarization. 

One may consider the possibility of using the temperature power spectrum to tighten 
constraints on t(z). The TT power spectrum is damped by a factor ~ e~ 2r ^ on sufficiently 
small scales due to free electrons scattering with microwave photons. The WMAP TT 
power spectrum is cosmic variance limited for I < 548 [5J and thus, the effect of dark matter 
annihilation is potentially observable in the TT spectrum. The top panel of Fig. [4] shows 
the TT power spectrum (red points) with no dark matter annihilation taken into account, 
compared with dark matter Model #3 (m x = 10 GeV, c = 15) shown by dashed lines. As 
expected, the power is slightly reduced when dark matter annihilation is included, due to 
the larger optical depth. However, this damping is nearly indistinguishable from the effect 
of reducing the overall normalization of the power spectrum set by A s , the amplitude of 
the primordial scalar power at the end of inflation. The bottom panel of Fig|4] compares 
the TT power spectrum with A s lowered by ~ 5% to the dark matter model, showing very 
good agreement. Reducing A s would result in slightly lower power on large scales in the 
TT spectrum (which are unaffected by reionization at high redshifts), but this difference is 
hindered by the larger cosmic variance at small /. 

The correlation between t(z) and A s may be quantified by means of the correlation 
parameter p constructed from the inverse of the Fisher matrix. For a 2 x 2 Fisher matrix, 
one readily finds: 

P = -^== (16) 



where 



I 



SQ 6Q 



5Ci represents the la error in C\ and the commas indicate a derivative w.r.t the parameter 
being varied. For an assumed form of C\ ~ A s e~ 2r<yZ \ we find p — 1, implying the Fisher 
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FIG. 4: The TT power spectrum. Partial ionization by dark matter results in a larger optical 
depth, causing a suppression of power (top panel). Decreasing the amplitude of the primordial 
scalar power A s by ~ 5% can mimic the effect of t(z)) (bottom panel). 

matrix is singular, with infinitely large errors in r(z) and A s . In practice, p is not exactly 1 
and the Fisher matrix is invertible. One might also ask whether alternate scenarios may be 
distinguishable from the effect of dark matter annihilation. To test this, we constructed a 
reionization history by modifying the recombination coefficient as follows: 



ttncw(^) = Oc{z) 



l + Z 

l + z 



(18) 



at redshifts below zq « 800, i.e. when the residual electrons start to become sub dominant. 



a(z) is the standard recombination coefficient (see Eq. (|6j)). Thus, Eq. (18) considers an 
ionization history with a faster fall off of the recombination coefficient with redshift. We 
calculate p by considering a simple 3 parameter Fisher matrix consisting of t(z), A s and 
spectral index n s , and marginalizing over n s . More parameters need to be included for a 
careful study, but this simple analysis suffices to show the degree of correlation between the 
parameters. We consider two cases with excess optical depth: (i) dark matter annihilation 
at high redshifts (with m x = 50 GeV, no halos) with the standard recombination coefficient, 



and (ii) no dark matter, but a smaller recombination coefficient as in Eq. (18) (with p 



0.985). For the TT power spectrum, we find for case (i), a correlation between t(z) and A s 
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FIG. 5: The large angle EE (top panel) and TE (bottom panel) power spectra. The la cosmic 
variance error bars are drawn for baryonic Model A. The solid (black) curve is drawn for the 
baryonic model with reduced A s . The dashed curve (blue) is plotted for dark matter Model #3 
(m x =10 GeV, c=15). It is harder to distinguish different models using the TE power spectrum 
because Cf E oc t(z), while C EE oc t(z) 2 . Following the WMAP convention, the bottom panel 
shows (I + l)Cf E /2-K and not 1(1 + l)Cf E '/2tt 



of p ~ 0.9989. For case (ii), we find p ~ 0.9993. The la errors in r(z) given by J[F 1 ]xi 
are found to be ~ 0.018 for case (i) and ~ 0.023 for case (ii). The exact values of p and 
Sr(z) will be different when more parameters are included. However, this simple calculation 
shows that r(z) and A s are strongly degenerate. We also see that the effects of dark matter 
annihilation cannot easily be distinguished from an alternate recombination theory. 

Fig. [5] shows the corresponding polarization power spectra. We verify that reducing A s 
has almost no effect on the large angle polarization power. Thus, while the dark matter 
model cannot be separated from the standard scenario using the TT power spectrum when 
A s is variable, they are easily separable using EE polarization data provided instrument 
noise is not a limiting factor. Also shown is the large angle TE power spectrum. One sees 
that the EE power spectrum is more useful than the TE power spectrum for distinguishing 
models, provided quality EE data exists for / < 20. The contribution due to dark matter 
halos mostly affects multipoles / ^ 20, while annihilation of free dark matter at high redshifts 
influences the multipoles / > 10. An ideal, cosmic variance limited observation would be 
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able to place constraints on both the halo model, as well as particle physics parameters. 
The Planck experiment is expected to be cosmic variance limited up to I ~ 10 [31] and have 
good signal to noise up to Z ~ 20 [30]. The CMBPol experiment is expected to perform even 
better. These future experiments may be able to distinguish between different reionization 
models, and place constraints on dark matter halos as well as particle parameters. 

IV. CONCLUSIONS. 

In this article, we have looked at how dark matter models influence the CMB polarization 
power spectrum, and whether current, and future observations can be used to distinguish 
dark matter reionization models from the standard theory. 

In Section II, we briefly reviewed dark matter annihilation resulting in ionization and 
heating, and derived an expression for the energy absorbed per gas atom £,(z). We also con- 
sidered inverse compton scattering of high energy charged particles with the CMB, resulting 
in lower energy photons (Fig. 1). When the inverse compton process is efficient, can 
be expressed in a much simpler form. At very high redshifts, £(z) can be further simplified 
since most of the released energy is absorbed close to the redshift of particle annihilation. In 
this regime, the ionization directly probes / cm (cr a t>) /m x , where / em is the fraction of energy 
transported by electromagnetic processes. In Section III, we discussed the optical depth 
contribution at different redshifts and the resulting CMB polarization power spectrum. We 
first considered dark matter annihilation at high redshifts z > 100, both with and without 
residual electrons (Fig. 2, Table 1). Taking residual electrons into account makes dark mat- 
ter annihilation less effective except for very small particle masses. We then considered dark 
matter annihilation at all redshifts, and calculated the ionization history for different dark 
matter masses, and concentration parameters (Fig. 3). We also considered the temperature 
power spectrum, and showed that it is not as useful as the polarization power spectrum in 
distinguishing different reionization scenarios (Figs 4,5). The EE power spectrum is more 
sensitive to reionization than the TE power spectrum. 

In this article, we have restricted ourselves to the case of a thermal relic with cross 
section (a a v ) = 3 x 10~ 26 cm 3 /s. However, dark matter annihilation at high redshifts probes 
fem(o'a.v)/m x and thus, our results are valid for a wide range of masses, for corresponding 
annihilation cross sections. We considered the effect of different particle masses, as well as 
different concentration parameters. The halo contribution which is dominant for z < 60 
affects the low multipoles / ^ 20 in the polarization power spectrum, while free dark matter 
annihilation at redshifts z > 60 primarily affects multipoles / > 10. We first looked at 
DM Model #1 with m x = 100 GeV and c = 15, and compared it with a 2-step standard 
reionization model (Model A). The effect of dark matter is not seen in the large angle 
polarization, even with cosmic variance limited data up to / = 30. Dark matter annihilation 
becomes important when the particle mass is decreased. We considered two models with 
m x = 10 GeV, with concentration parameters c = 5 and c = 15 (DM Model #2, and DM 
Model #3 respectively). With cosmic variance limited data up to / ~ 20, DM Model #2 
and DM Model #3 may be distinguished at high significance from Model A. Thus future 
CMB polarization measurements might allow us to put constraints on light (on the weak 
energy scale) dark matter candidates and on models with an enhanced annihilation cross 
section. Current large angle CMB polarization measurements are consistent with a standard 
baryonic reionization theory. Future experiments such as Planck would be able to shed more 
light on the physics of reionization. 



14 



Acknowledgments 

We thank Nishikanta Khandai, Roberto Trotta, Licia Verde, and Eiichiro Komatsu for 
helpful discussions. A.N. thanks the Bruce and Astrid McWilliams Center for Cosmology 
for financial support. D.S. was supported by the Deutsche Forschungsgemeinschaft under 
grant GRK 881. 



[1] A. Kogut et al., for the WMAP collaboration, Astrophysical Journal Supplement Series, 148, 
161 (2003) 

[2] L. Page et al., for the WMAP collaboration, Astrophysical Journal Supplement Series, 170, 
335 (2007) 

[3] M.R. Nolta et al., for the WMAP collaboration, Astrophysical Journal Supplement Series, 
180, 296 (2009) 

[4] J. Dunkley et al., for the WMAP collaboration, Astrophysical Journal Supplement Series, 
180, 306 (2009) 

[5] D. Larson et al., for the WMAP collaboration, |arXiv:1001.4635| ^l [astro-ph.CO] (2010) 

[6] D. Spolyar, K. Freese, P. Gondolo, Phys. Rev. Lett., 100, 051101 (2008); F. Iocco, Astrophys. 

J., 677, LI (2008); K. Freese, D. Spolyar, A. Aguirre, JCAP, 11, 14 (2008); K. Freese, P. 

Bodenheimer, D. Spolyar, P. Gondolo, Astrophys. J. 685, L101 (2008); F. Iocco, A. Bressan, 

E. Ripamonti, R. Schneider, A. Ferrara, P. Marigo, Mon. Not. R. Astron. Soc, 390, 1655 

(2008); A. Natarajan, J.C. Tan, BW. O'Shea, Astrophys. J. 692, 574 (2009); D. Spolyar, P. 

Bodenheimer, K. Freese, P. Gondolo, Astrophys. J., 705, 1031 (2009) 
[7] E. Pierpaoli, Phys. Rev. Lett. 92, 031301 (2004), X. Chen, M. Kamionkowski, Phys. Rev. D 

70, 043502 (2004), L. Zhang, X. Chen, Y.A.Lei, Z.G.Si, Phys. Rev. D 74, 103519 (2006), M. 

Mapelli, A. Ferrara, E. Pierpaoli, Mon. Not. R. Astron. Soc. 369, 1719 (2006), E. Ripamonti, 

M. Mapelli, A. Ferrara, Mon. Not. R. Astron. Soc. 375, 1399 (2007); 374, 1067 (2007), M. 

Mapelli and E. Ripamonti, Mem. Soc. Astron. Ital. 78, 800 (2007), L. Chuzhoy, Astrophys. 

J. 679, L65 (2008) 
[8] N. Padmanabhan, D.P. Finkbeiner, Phys. Rev. D 72, 023508 (2005) 
[9] A. Natarajan, D.J. Schwarz, Phys. Rev. D 80, 043529 (2009) 
[10] A. Natarajan, D.J. Schwarz, Phys. Rev. D 78, 103524 (2008) 
[11] A.V. Belikov, D. Hooper, Phys. Rev. D 80, 035007 (2009) 
[12] J. Chang et al, Nature 456, 362 (2008). 

[13] O. Adriani et al. [PAMELA Collaboration], Nature 458, 607 (2009) |arXiv:08 10.4995] [astro- 
ph]]. 

[14] M. Cirelli, F. Iocco, and P. Panci, J. Cosm. Astropart. Phys., Issue 10, 009 (2009), G. Hiitsi, 
A. Hektor, M. Raidal, Astron. and Astrophys., 505, 999 (2009), Q. Yuan, B. Yue, X. Bi, X. 
Chen, X. Zhang, |arXiv:0912.2504"l [astro-ph] (2009) 

[15] T. Kanzaki, M. Kawasaki, K. Nakayama, |arXiv:0907.3985] [astro-ph] (2009); F. Iocco 
|arXiv:0912.1630| [astro-ph] (2009) 



[16] J. Kim, P. Naselsky, |arXiv:1001.4691| rl [astro-ph.CO] (2010) 
[17] Legacy Archive for Microwave Data Analysis, 

\protect\vrule widthOpt\protect\href {http : //lambda. gsfc .nasa. gov/product/map/dr3/parameter: 
[18] E. Komatsu et al., for the WMAP collaboration, arXiv:1001.4538> l [astro-ph.CO] (2010) 



15 



A. Lewis, A. Challinor, A. Lasenby, Astrophys. J. 538, 473 (2000) 

A. M. Green, S. Hofmann, and D. J. Schwarz, Mon. Not. R. Astron. Soc. 353, L23 (2004) 



A. M. Green, S. Hofmann and D. J. Schwarz, JCAP 0508, 003 (2005) |arXiv:astro- 
|ph/0503387| . 

J. M. Shull and M. E. Van Steenberg, Astrophys. J. 298, 268 (1985) 

T. Kanzaki, M. Kawasaki, Phys. Rev. D, 78, 103004 (2008); 

S. Furlanetto, S.J. Stoever |a?Xiv:0910.4410[ astro-ph] (2009) 

M. Valdes, C. Evoli, A. Ferrara, |arXiv:0911.Tl25 [astro-ph] (2009) 

R. Weymann, Phys. Fluids 8, 2112 (1965); 

S. Seager, D.D. Sasselov, D. Scott, Astrophys. J. 523, LI (1999), Astrophys. J. S. 128, 407 
(2000) 

F. C. Jones, Phys. Rev. 167 (1968); G.R. Blumenthal, R.J. Gould, Rev. Mod. Phys. 42, 237 
(1970) 

R. Schlickeiser, "Cosmic Ray Astrophysics", Springer- Verlag, Germany (2002); O. Petruk, 

Astron. Astrophys. 499, 643 (2009) 

The Planck Blue Book, ESA-SCI(2005)1 

M. Zaldarriaga et al, arXiv:0811.3918| [astro-ph] (2008) 

J.M. Shull, A. Venkatesan, Astroph. J. 685, 1 (2008) 

N.Y. Gnedin and J. P. Ostriker, Astrophys. J., 486, 581 (1997); G.P Holder, Z. Haiman, M. 
Kaplinghat, L. Knox, Astroph. J. 595, 13 (2003) 

J. F. Navarro, C. S. Frenk, and S. D. M. White, Astrophys. J. 490, 493 (1997) 

R.K. Sheth, G. Tormen, Mon. Not. R. Astron. Soc, 308, 119 (1999); R.K. Sheth, H.J. Mo, 

G. Tormen, Mon. Not. R. Astron. Soc, 323, 1 (2001) 



16 



